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ABSTRACT 

Although the highest redshift QSOs (z > 6.1) are embedded in a significantly neutral background 
universe (mass-averaged neutral hydrogen fraction > 1%) as suggested by the Gunn-Peterson absorp- 
tion troughs in their spectra, the intergalactic medium in their vicinity is highly ionized. The highly 
ionized region is generally idealized as spherical and called the Stromgren sphere. In this paper, by 
combining the expected evolution of the Stromgren sphere with the rule that the speed of light is fi- 
nite, we illustrate the apparent shape of the ionization fronts around the highest redshift QSOs and its 
evolution, which depends on the age, luminosity evolution, and environment of the QSO (e.g., the hy- 
drogen reionization history). The apparent shape may systematically deviate from a spherical shape, 
unless the QSO age is significantly long compared to the hydrogen recombination process within the 
ionization front and the QSO luminosity evolution is significantly slow. Effects of anisotropy of QSO 
emission are also discussed. The apparent shape of the "Stromgren sphere" may be directly mapped 
by transmitted spectra of background sources behind or inside the ionized regions or by surveys of 
the hyperfine transition (21 cm) line emission of neutral hydrogen. 

Subject headings: cosmology: theory — early universe — galaxies: high redshift — quasars: general 
— quasars: absorption lines — relativity 



1. INTRODUCTION 

The existence of the Gunn-Peterson absorption trough 
in the spectra of the highest redshift QSOs (z > 6.1) sug- 
gests that hydrogen in the early uni verse is significa ntly 
neut ral (> 1% in mas s avera ge; e.g.. | Fan et alJl2 002t see 
also iWvjthe fc Loebl l2004al iMesineer fc Haimanl 120041 
lYu fe LuTl2005[r However, hydrogen in the vicinity of 
the QSOs may be highly ionized by ioni zing photons 
from QSOs (i.e., the proximity effect; e.g.. iBaitlik et al.1 
1988), as indicated in the QSO spectra by the transmis- 
sion of flux at wavelengths between the Gunn-Peterson 
trough and the Lya line center. These highly ionized re- 
gions around the QSOs are generally idealized as spher- 
ical in the rest frame of the QSO and called Stromgren 
spheres. 2 The expansion of the Stromgren sphere de- 
pends on the QS O luminosity evolution and its envi- 
ronm ent (e.g.. lDonahue fc ShulllllQS^: iShapiro fc Girouxl 
I19S7T) . Observers on Earth perceive an object (here the 
highly ionized region) by the light reaching Earth at the 
same moment. Since the speed of light is finite, the ap- 
parent shape of the expanding ionization front around 
t he QSO may not b e spher ical. 

iWvithe fc Loebl l)2004b|) show that the ionization 
fronts along the line of sight and transverse to it may 
have different observed sizes ( see also discussions on the 
size along the line of sight inlWhite et al.ll2003l and re- 
lated discussions in lPentericci et al.ll2002|) . In this paper 
we generalize the time-delay effect to all angles around 
the highest redshift QSOs and illustrate the expected ap- 
parent shape of the surrounding ionization front and its 

1 Hubble Fellow 

2 These highly ionized HII regions (or Stromgren spheres) are 
well defined only for QSOs that have the Gunn-Peterson trough. 
QSOs with relatively low redshifts do not have an unambigu- 
ous Stromgren sphere because their environment has already been 
highly ionized by stars (e.g., sec Adclbcrgcr 2004). 



evolution (for the relativistic time- delay effects in other 
astrophysical contexts, e.g., see also ICouderHlT989t iReTsI 
1966). We show how the shape depends on the QSO 
luminosity evolution and its environment (such as the 
neutral hydrogen fraction and dumpiness). Comparison 
of observations with the model prediction of the appar- 
ent shape may provide constraints on QSO properties, 
the reionization history of the universe, and/or the cos- 
mo logical parameter 17a ( by the Alcock-Paczyhski test; 
lAlcock fc Paczvhskilll979fl . 

This paper is organized as follows. In § [2 we illustrate 
the geometry of the time-delay effect on the apparent 
shape of the ionization front around a QSO. In §0 we 
study the evolution of the Stromgren sphere in the QSO 
rest frame. We combine the evolution and the time-delay 
effect to illustrate the evolution of the apparent shape of 
the ionization front in § 0] The effects of anisotropy of 
QSO emission are discussed in §0 Possible observational 
tests of the apparent shape are discussed in § El Our 
conclusions are summarized in § 

In reality the boundary of an ionized region may not 
be exactly round even without considering the relativis- 
tic time-delay effect and the QSO anisotropic emission. 
Since the Stromgren sphere is idealized as spherical by ig- 
noring other possible random fluctuations on the bound- 
ary, the apparent shapes shown in this paper are also 
idealized as smooth and may more appropriately repre- 
sent the average over the fluctuations. 

2. GEOMETRY 

We use Figure n to illustrate the geometry of the ap- 
parent shape of the ionization front around a QSO at 
point O. For simplicity, we ignore the Hubble expan- 
sion within the ionization front, since here the physi- 
cal size of the ionized region around the highest redshift 
QSOs is much smaller than the Hubble scale and the 
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characteristic timescale of the QSO luminosity evolution 
is much smaller than the age of the universe at z ~ 6. 
The derivation below is a simple application of the rule 
that the speed of light is finite. Suppose that the nu- 
clear activity of the QSO is triggered at cosmic time ti- 
The ionizing photons emitted from the QSO will then 
ionize surrounding neutral hydrogen. For simplicity, we 
assume that the QSO emission is isotropic here and dis- 
cuss the anisotropic effect of the emission in § In the 
rest frame of the QSO, a spherical ionization front will 
form and expand (see the dotted circles in Fig. , sep- 
arating the inside highly ionized region and the outside 
neutral or partly neutral region. We denote the physical 
(proper) radius of the expanding sphere at cosmic time 
t (> ti) by R(t;ti). The detailed evolution of R(t) is 
related to the QSO luminosity evolution and its environ- 
ment, which will be determined in § [21 The observer on 
Earth perceives the ionization front by the light reach- 
ing Earth at the same moment. If photons emitted from 
(or passing through) the ionization front, e.g., point A 
in Figure reach the distant observer at the same time 
as photons from the QSO, the photons from A should be 
emitted from (or pass through) A at cosmic time 

t A =t(z Q ) + R(t A ) cos 9/c, < 9 < 7T, (1) 

where zq is the QSO redshift, t(zQ) is the cosmic time 
corresponding to zq, c is the speed of light, and 9 is the 

angle between OA and the observer's line of sig ht od 
(see Fig. 0. For QSO photons to reach the ionization 
front A at cosmic time t A , they should emit from the 
QSO O at an earlier time 



to=t A - R(t A )/c. 



(2) 



Given any to (U < to < t A ), the solution of equa- 
tion J21 t A is unique since the expansion speed of R{t A ) 
is slower than c. Thus we may define such a function 
r(r; ti) = R(t A ; ti) by changing the independent variable 
from t A to t, where r = to — ti is the age of the QSO at 
to- According to the above, the apparent shape of the 
ionization front (r, 9) should satisfy the following equa- 
tion: 



•(t)(1 



cos 



9)=c(r Q -r), 



< 9 < 



(3) 



where tq = t(zQ) — ti is the age of the QSO at cos- 
mic time t(zQ). The apparent shape (illustrated by the 
solid curve in Fig. is rotationally symmetric along the 
observer's line of sight to the QSO. The shape is gener- 
ally not a sphere [unless r(r) does not change with r]. 
According to equation 10, given tq and an increasing 
function r(r), the apparent size r decreases with increas- 
ing 9. Because of the existence of the term "— r" in the 
right-hand side of equation |J3}, the apparent shape (r, 9) 
is neither a spheroid (the apparent shape of relativisti- 
ca lly expanding sph erical shell of radio sources discussed 
in lReesl[l966l Il967j) unless r(r) is a linear function of r 
nor a "light-ech o" paraboloid around a nova or supernova 
(|Coudercl ll939TK 

We define the apparent radial expansion speed of the 
ionization front at a given 9 by u app = dr' /dt', where 
r' = (1 + z)r is the comoving distance of r and dt' = 
(1 + z)drQ is the time interval measured by the observer 
with including the effect of the cosmological time dila- 
tion. According to equations 01-10, the apparent radial 



expansion speed of the ionization front in units of c at a 
given 9 is given by 

dr/dr 



^app 



app 



1 dr 
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dR/dt A 
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(4) 



(5) 



c — (dR/dt A ) cos 9' 
Note that the physical expansion speed of the ionization 
front in the QSO rest frame is 

dR/dt A = (dr/dr) /(l + cr 1 dr/dr) < c, (6) 
which is consistent with special relativity that real infor- 
mation does not travel faster than the speed of light. If 
dR/dt A is highly relativistic, the apparent expansion is 
super luminal at sufficiently small 9, i.e., /3 app > 1, and 
/3 app is close to 1/2 when 9 — > it. We can also define 
the apparent transverse expansion speed in units of c 
through /3 app ,_L = /?appSin6>, which gives the convention- 
ally transverse superluminal phenomenon at sufficiently 
small 9 if dR/dt A — > c. According to equation 10, the 
lower limit on the QSO lifetime can be constrained sim- 
ply by the apparent sizes of the ionization front, e.g., 
tq > 2r(9 = n)/c or r Q > r(9 = tt/2)/c. 

3. THE EVOLUTION OF r(r) 

We denote the mean neutral hydrogen fraction in the 
region surrounding the QSO prior to the triggering of 
its nuclear activity by Xm = (nm) I (^h) (~ several per- 
cent to 10-20% in lYu fc Lul 120051: and the rest of the 
hydrogen has been ionized b y stars or co Uisionall y ion- 
ized in large halos; als o cf.. IMesinger fc Haimanl 120041 
or for a higher fraction IWvithe 7TLoebTl2004a|) . where 
(riHi) is the average proper number density of neutral 
hydrogen and (nij) is the density of total hydrogen. Af- 
ter the nuclear activity is triggered, the expansion of the 
ionization front (or the Stromgren sphere) in the rest 
frame of the QSO as a function of r (not as a function 
of t A here) may be desc ribed by the following eq uation 
Ce.g.. iDonahue fc Shul| Il987t iShapiro fc Gircm3 1198 ^ 
iCen fc Haimanl l2000t llVIadau fc Reesl 120001: lYTfcL 
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3 dr ' 3 

where N-phs(r) is the ionizing photon emission rate of the 
QSO at its age r, a B (= 2.6 x 10~ 13 cm 3 s" 1 at T = 10 4 K) 
is the hydrogen recombination coefficient to excited lev- 
els of hydrogen, C = (^Hi)/( n n) 2 is the clumping factor 
describing the effective dumpiness of hydrogen ionized 
by QSO photons, and "(■ • •)" represents the space av- 
erage over the sphere with radius r. Note that here C 
is a little different from the conventional definition of 
(ng)/(riH) 2 - The second term on the right-hand side of 
equation 10) accounts for the recombination of ionized 
hydrogen within the Stromgren sphere. In equation (0 
we set the parameters (xhii (^h), C) to be the values at 
cosmic time t = t, t + r, and their changes during the 
time interval r/c of the QSO photon propagation in the 
Stromgren sphere have been assumed to be negligible. 
For simplicity, we also neglect their dependence on r and 
set them to be constants in the calculations below. Thus, 
the solution of equation 10 is given by 

1/3 
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Fig. 1. — Schematic diagram for the apparent shape of the ion- 
ization front around a QSO (with the Gunn-Peterson trough in 
the spectrum). The QSO is at point O and the distant observer 
on Earth is located in the direction of OC. The QSO emission is 
assumed to be isotropic. The dotted circles represent the ioniza- 
tion fronts at a given cosmic time (t a and t g for inside and outside 
circles, respectively) in the rest frame of the QSO. Hydrogen in- 
side the circle (or the Stromgren sphere) is highly ionized by QSO 
photons. The solid curve is the apparent ionization front, which 
follows equation Photons emitted from (or passing through) 
the apparent ionization front (e.g., from point A at tA and B at 
tg shown in this figure) reach the observer at the same time. The 
apparent shape of the ionization front is rotationally symmetric 
along OC and is generally not a sphere/spheroid/paraboloid. The 
apparent shape could be mapped by observational spectra of back- 
ground objects (e.g., S in the figure) or the tomograph of HI 21 cm 
line emission. 



where the recombination timescale 
Tree = xm (C (n H ) cub ) ~ 1 
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l + z 



(9) 



characterizes the role of the recombination term in equa- 
tion J7J). 

The physical expansion of the ionization front in 
the rest frame of the QSO dR/dtA was discussed in 
I White et alJ '|2( )()•'>• by accounting for finite light travel 
time and ignoring the rec ombination terms , and w as gen- 
eralized in equation (1) of lWvithe fe Loebl lj2004aj) . Here 
applying equation (JJJ into equation ijfJJl, we obtain the 
physical expansion of the ionization front as follows: 



dR 



iV phs - f a B C(n a ) 2 R 3 



dt A 4TTcxm{nii)R 2 + N phs 



(10) 



we have r = tq (eq.|3J) and thus N p hs is simply N^^tq); 
at 9 = it/2 we have r = tq — R/c and thus iV p h s is simply 

Nph s (TQ — R/c). Note that in principle the apparent ion- 
ization front at 9 = tt/2 may not necessarily correspond 
to the maximum angle of the apparent ionized region 
extended on the sky (see Fig. 

Note that in equation J7J N p hs(T~) is a non-decreasing 
function of r so that r(r) expands with increasing r. 
When iVphs(T) decreases with increasing r (e.g., close to 
the quenching of the QSO nuclear activity) and r 3> r rec , 
the ionization front shrinks, and its evolution may not be 
described by equation @ if the characteristic decreasing 
timescale of iVphs(T) is significantly short compared to 
the recombination timescale. For simplicity, in this pa- 
per we ignore the decreasing phase of the QSO luminosity 
evolution, which has been shown to at least not dominate 
the main population of luminous QS Os observed at rel- 
atively low redshifts l)Yu fc Lull20"04|) . 

Below we assume two simple models of N p hs(i~) (see 
also lYu fc Lull2005|) : in model (i) iVp hs (r) = N phs>i is a 

constant, and in model (ii) Np^r) — iV p h s ,i ex P( T / T s) 
increases exponentially with r, where t$ — 4.5 x 
10 7 yr[0.1e/(l — e)] is the Salpeter timescale and e is the 
mass-to-energy conversion efficiency. When r <C r rcc in 
model (i) or r <C min{r rcc , rg} in model (ii), the effects of 
recombination and the evolution of the QSO luminosity 
are both negligible and r(r) is given by 
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Given equation (|1 1|) . the changing rate of r(r) with r is 
larger than c if r < [rs/{3cT rcc )] 3 / 2 T rcc (i.e., r(r) > 3ct 
here). As the ionization front expands, if r 3> T rcc (and 
Tree *C ts in model ii), the recombination within the 
highly ionized HII region is approximately balanced by 
the emission of the ionizing photons from the QSO, and 
we have 



r(r) 



+T S 



r/rs 



1/3 



for model (i), 
for model (ii), 



.(13) 



where iV phs = ^ P h s (r) = N phs {t A - U - R/c). At 6 = 



which is insensitive to the detailed value of r for model 
(i) and depends on r only through A^ p i 1s (t) in model (ii). 
In case that r rcc 3> ts in model (ii), r(r) at t > tj 
still follows equation l|13|l , but the effect of recombination 
is not important and r(r) is mainly determined by the 
QSO emission. Given equation (|13|l . for the exponential 
increase of the luminosity evolution in model (ii), the 
changing rate of r(r) with r is larger than c if r(r) > 
3cr s . 

4. THE EVOLUTION OF THE APPARENT SHAPE OF THE 
IONIZATION FRONT 

By combining equations |J3J and ©, we obtain the 
evolution of the apparent shape of the ionization front 
(r, 6) for different models of the QSO luminosity evolu- 
tion, which is shown in Figure [5] The evolution of the 
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Fig. 2. — Expected evolution of the apparent shape of the ion- 
ization front around a QSO. The geometric configuration of the 
QSO (located at the crossing point of the dotted lines) and the 
observer is the same as in Figure^ The QSO emission is assumed 
to be isotropic. Both of the axes are in units of cr rec (see T rec 
in eq. Top panels (a)-(b) are for model (i) of TVp^s (see S IH. 
and middle and bottom panels for model (ii) with Ts/r Tec = 10 
(c-d) and 3 (e-f), respectively. The parameters (JV p h s ,ii ( n H)i C) 
are set so that r(r = 10r rcc )/(cr rcc ) = 4 in the left panels (a,c,e) 
and 8 in the right panels (b,d,f). In each panel, different curves 
represent the apparent shapes with different tq / r ICC ■ The tq / r Tec 
increases from inward to outward with 5 log(rQ /t tcc ) = 0.25, and 
the outmost curve has tq/t icc = 10. 



apparent shape and size depends on the QSO luminosity 
evolution and its environment (e.g., the hydrogen reion- 
ization history or structure formation) through the pa- 
rameters (tq, ts, T r ec, rg). In Figure [SJ the parameters 
(N p h s ,i,Xm,C, (nn)) are set so that r s /{cT lcc ) is larger 
in the right panels than in the left ones. Top panels (a) 
and (b) are for model (i) of iV p h s (■?")• As seen from them, 
the apparent shapes of the ionization fronts are generally 
not a sphere, especially if tq is short (inner curves) and 
rs/{cT ICC ) is large (e.g., because N phSyi is high or x m is 
low; right panel). The apparent shape expands as tq 
increases. When tq is long enough, the apparent shape 
at small 9 (see the definition of 9 in Fig. ^) starts to 
be almost independent of tq and approaches a sphere, 
since the radius of the Stromgren sphere approaches the 
equilibrium value rs (eq. 113(1 at which the ionizing pho- 
ton emission is balanced by the recombination. With in- 
creasing tq , the shape at large 9 continues to expand and 
finally the whole apparent shape becomes spherical (see 
the outmost curve in panel a). According to equation 
<|12|) . the stronger the iV p h s ,i or the smaller the C, the 
larger the final equilibrium sphere radius. In addition, 
in case that rs/(cT rcc ) is significantly small (<C 1), the 
expansion may become non-relativistic at a time much 
earlier than tq ~ r rec , thus the apparent shape may ap- 
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Fig. 3. — Expected apparent radial expansion speed of the ion- 
ization front in units of c as a function of tq / T roc (eq. [IJ . Panel 

(a) (model i of N p hs) na s the same parameters as those in Fig- 
ure Ha. and panel (b) (model ii) the same as those in Figure ITI:. 
In each panel, different curves give the speeds at different 8, which 
increases from to 180° with 88 = 30° from top to bottom at the 
low- tq /r r ec end. The apparent speed at small 8 is super luminal at 
small tq I T roc in both panels and also at large tq / T rec for model 
(ii) in panel (b) . For simplicity, we do not show the apparent trans- 
verse expansion speed /3 apPi ^ = /3 a ppsin0 here, which also gives 
the similar superluminal phenomenon except at 8 = 0. 



proach a sphere even if tq -C r rac (but tq > t(tq)/(3c)). 
Figure [5]>f illustrate the apparent shapes for model (ii) 
of iVphs(T), where the timescales of increasing luminos- 
ity are set to have rs/r rcc = 10 and 3, respectively. In 
model (ii) the size of the apparent shape is sensitive to 
tq even if tq is long because iV p h s increases with time 
exponentially; and the shape does not evolve towards a 
spherical shape, especially if t$ is small (bottom panels 
e-f). 

For the realistic case of the highest redshift QSOs with 
Gunn-Peterson troughs obse r ved so far (z ~ 6.1 — 6.4; 
iF^EfflllWhite et al.ll200l . IYu fc Lil ll20f)l study the 
environmental effect of their surrounding regions and 
obtain xhi ~ 0.1 and C ~ 30-40 (which are shown 
to be consistent with the observed Stromgren radii). 
Then the recombination timescale r roc is ~ 4 x 10 6 yr 
(eq. IHJ , which is signific antly shorter th an the QSO life- 
time (> 4 x 10 7 yr; e.g. JYu fc Lull2004|) . Thus the ages 
of most observed QSOs should be long enough that the 
sizes of the ionization fronts along the line of sight are 
close to the equilibrium values. For the observed QSOs 
with N p yi S {tq) ~ 10 57 s _1 , t he apparent siz e at 9 = is 
t-(tq) - 5Mpc(- 3-4cT roc ) l|Yu fc Lul2005j) . and the ap- 
parent shapes should be close to the outer curves shown 
in Figure and c (see also Fig. 0] if the QSO emission 
is anisotropic). Statistically a small fraction of observed 
QSOs would also have significantly short ages and the ap- 
parent shapes of their ionization fronts should systemati- 
cally devi ate from a spherical shap e (see the inner curves 
in Fig. I2J). W vithe fc Loebl l(2004a|) obtain a much higher 
Xhi and a low QSO age from the observed Stromgren 
radius of a QSO by using a much smaller clumping fac- 
tor; in this case, the recombination timescale is signifi- 
cantly long (up to 10 9 yr) and t(tq) is not significantly 
smaller than 3ctq, so the apparent shapes of most ob- 
served QSOs should be more likely to systematically de- 
viate from a spherical shape. 

Using equation JSJ, we illustrate the apparent expan- 
sion speed of the ionization front as a function of tq in 

Figure|ni Panel (a) is for model (i) of iV p h s and panel (b) 
for model (ii). As seen from Figure |3 the speed is su- 
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per luminal at small 9 and tq / T rec . For the specific case 
of 6 = 0, we have t = tq and /3 app = c~ x dr /dr (eqs. |3] 
and , and the corresponding superluminal expansion 
is also mentioned or discussed in some recent literature 
llOen fc Ha,ima,nl2nnflHWhite et a,l.l2003HWvithe fc Loehl 
2004b). In panel (b), the speed at small 9 may also 
become superluminal even at large tq/t vcc because the 
ionizing photon emission rate and r(r) (ea. H3fl increases 
exponentially in model (ii) . The expansion of the highly 
ionized regions around the highest rcdshift QSOs with 
Gunn-Peterson troughs provides another astrophysical 
example of the superluminal phenomenon (for other ex- 
ampl es, e.g., see reviews in Blandford. McKee fc Reesl 
Il977t iFalla. Flovd. fc Potterll2003j) . 

5. EFFECTS OF ANISOTROPY OF QSO EMISSION 

In the above, QSO ionizing photons have been assumed 
to be emitted isotropically. However, in reality the QSO 
emission may be anisotropic at least in these two cases: 
(i) a torus is surrounding the central black hole in the 
QSO, which blocks the ionizing photons from the central 
engine; (ii) the UV ionizing photon emission from the 
accretion disk surrounding the black hole is anisotropic 
(see also di scussions o f th e effects of anisotropy o f QSO 
emission in lOenll2003l and iWvithe fc Loebll2004bfl . 

In case (i), the highly ionized region in the rest frame 
of the QSO may appear simply like a large-scale version 
of the [OIII] ionization "cone" detected in local Seyfert 
2 galaxies. The expected apparent shape is easily corre- 
spondingly adjusted depending on the normal direction 
and the opening angle of the torus. 

In case (ii), we assume that the emergent intensity from 
the accretion disk is oc (1 + 2\u\) because of the limb- 
darkening effect (e.g. . lKrolikll999|) . where fi is the cosine 
of the angle between the disk normal and the emitting di- 
rection (so far there is no much evidence to show that UV 
photons from QSOs are strongly beamed). Here the QSO 
anisotropic emission is symmetric about the disk plane 
and axisymmctric about the disk normal. Although both 
the relativistic time-delay and the anisotropy of the QSO 
emission may make the apparent shape deviate from a 
spherical shape, their effects arc different. As mentioned 
above, the time-delay effect decreases r(r) with increas- 
ing 9 (see eq. |3J) and does not break the rotational sym- 
metry of the ionized region about the observer's line of 
sig ht od. For the anisotropic effect here, the detailed 
dependence of the apparent size on 9 is significantly af- 
fected by the detailed geometric configuration between 
the accretion disk and the observer. The rotational sym- 
metry of the apparent shape around the observer's line 
of sight is generally broken unless the disk normal is co- 
incident with the observer's line of sight. We denote the 
angle between the disk normal and the observer's line of 
sight by 9' (0 < 6' < it/ 2). For 9' jt 0, the cross sections 
of the apparent ionized region are illustrated in Figure^] 
where the observer's line of sight is on the cross section, 
4> (0 < 4> < 7r ) is the angle of the cross section anti- 
clockwise from the plane in which ud and the disk nor- 
mal are located, and = sin 9 cos sin 9' ± cos 9 cos 9' \ . 
Although in general r(r) does not necessarily decrease 
with increasing 9 if the QSO emission is anisotropic, the 
apparent size in any direction from the central QSO with 
9 < 7r/2 is still systematically not larger than the size in 




Fig. 4. — Cross sections of the apparent ionized region af- 
ter including the anisotropy of the accretion disk emission due to 
the limb-darkening effect in Figure Ek . The solid lines with ar- 
rows in (a) and (b) are the disk normal passing the QSO. The 
emission is assumed to be symmetric about the disk plane and 
rotationally symmetric about the disk normal. The parameters 
(^Vphs,ii ("h)^) are set so that r[T(10r Tec , 9 = 0] = 4. After in- 
cluding the anisotropic effect, the rotational symmetry about the 
observer's line of sight mentioned in Figure^is broken, so the cross 
sections arc different in the top, middle, and bottom panels. As in 
Figure for each curve, the apparent size in any direction from 
the central QSO with 9 < tt/2 is not smaller than the size in the 
opposite direction due to the time-delay effect. See details in S 151 



the opposite direction because of the time-delay effect 
and the symmetry assumed for the anisotropic emission 
here. Given the detailed discussions in § if the QSO 
age is significantly long compared to the hydrogen re- 
combination process and the QSO luminosity evolution 
is significantly slow, the apparent ionized region follows 
the same symmetry as the QSO emission (i.e., symmetric 
about the disk plane and axisymmetric about the disk 
normal here). As seen from Figure the smaller the 
TQ/r rec and the smaller the 9', the less important the 
anisotropic effect. 

6. POSSIBLE OBSERVATIONAL TESTS 

We use Figurento illustrate the possible observational 
tests of the apparent shape. The discussion below may 
be generalized to the anisotropic case of QSO emission. 

If there exists a background object S behind the 
highly ionized region (see Fig. significant transmis- 
sion of Lya (and/or Ly/3) photons may be present on 
the Gunn-Peterson trough in the observational spectrum 
of S at wavelengths corresponding to redshifts z{tA) 
to z{ts) [somewhat like (but it may be much broader 
than) the observed emission spike on top of the Gunn- 
Peterson trough o f SDSS J1148+5251; see Figs. 4 or 5 in 
I White et al.ll200 3] . Similarly, the spectra of background 
objects lying inside the highly ionized re gions (see also 
discussions for inside objects in lCenl2003[) may also show 
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transmission features of Lya photons at the wavelength 
corresponding to the redshift when its light travels out 
of the ionized region. If the background object is the 
QSO itself, we have r = tq at 9 = (eq. |3J) and the size 
of the apparent ionization front \OC\ = t(tq), which 
can be measured directly from the wavelength difference 
between the Gunn-Peterson tro ugh and the Lya li ne cen- 
ter in the QSO spectrum fe.g. JWhite et al.ll2003|) . Note 
that the existence of Lya emitting galaxies lying behind 
or inside the highly ionized regions of the highest redshift 
QSOs is likely and detection of them is possible by cur- 
rent techniques, since the apparent "Stromgren sphere" 
with a physical size ~ 5 Mpc at z ~ 6 covers an area 
about 14 x 14 arcmin 2 on the sky and is large enough for 
this compared to the survey area s of some Lya surveys 
(see Table 1 in IStern et alJl2004]) that have successfully 
discovered galaxies at z ~ 6.5. By obtaining observa- 
tional spectra of a number of background sources with 
different projecting distances (R± = \OA\sin8) to the 
highest redshift QSO(s), one can measure the proper dis- 
tance of \AB\ (or \OA\, \OB\) for each source behind the 
region(s) (or \OB\ for inside sources) and then statisti- 
cally map the three-dimensional apparent shape of the 
highly ionized region(s) surrounding the QSO(s). The 
distribution of these distances as a function of R± may 
be used to statistically constrain the evolution of the ap- 
parent ionization fronts (and hence the QSO luminosity 
evolution and its environment), as illustrated in Figure[5J 

In addition, the apparent shape of the highly ion- 
ized region surrounding the QSOs with Gunn-Peterson 
troughs might also be revealed by radio surveys of the 
hyperfine transition 21 cm line emission of neutral hydro- 
gen (e.g., the Low Frequency Array, the Square Kilome- 
ter Array, the PrimevAl Structure Telescope, the Mileura 
Wide-Field Array) , a possible technique to map the neu- 
tral hy drogen distribution at h igh redshift (see discus- 
sions in lWvithe fc Loebll2r)04bT) . 

For the observer on Earth, the physical distance of R± 
can be obtained by R j_ = DaS4>, where Da is the angular 
diameter distance to the QSO and 5(f) is the angle of the 
QSO and the ionization front (e.g., point A) separated 
on the sky. The transformation of the angular and the 
redshift space to real space (R± and \AB\, \OA\ or \OB\ 
etc.) depends on the cosmological parameters. If the 
QSO age and environment are effectively constrained by 
the apparent ionization front (through \OC\, \AB\/\OC\ 
etc.) at small R± ~ and if the anisotropy of QSO emis- 
sion is effectively constrained by some cross sections of 
the apparent ionized region, the apparent physical sizes 
of the ionized region at other directions may be predicted 
by using the method in this paper. Thus, the ratio of the 
redshift difference between the QSO and the apparent 
ionization front at R± ~ to the maximum angle(s) of 
the apparent ionized region (or cross sections) extended 
on the sky might be used to provide constraints on the 
cosmological parameter Oa by the Alcock-Paczyhski test 
l)Alcock fc Paczvhskilll97ft) . since the angular size signif- 
icantly depends on Q\ but the redshift size does not at 
z ~ 6. 

7. CONCLUSIONS 

In this paper we have illustrated the apparent shape of 
the ionization front around the highest redshift QSOs, 
and its evolution, with Gunn-Peterson troughs in the 




Rj/IOCI Rj/|OC| 

Fig. 5. — Expected ratio of |AB|/|OC| of the apparent ionization 
front as a function of R±/\OC\ (see Fig.^for the definition of the 
proper distances of |AB|, |OC|, and R±). The parameters are 
the same as those in Figure [5] This figure illustrates that the 
evolution of the apparent ionization fronts around QSOs may be 
inferred from the distribution of |AB|. For observational methods 
of determining \AB\ and |OC|, see discussion in § |H| A similar 
study may also be done for \OA\ and |OB|. 



spectra. We have shown that the apparent shape and its 
evolution depend on the age and the luminosity evolu- 
tion of the QSO, the recombination process of ionized hy- 
drogen surrounding the QSO, and the anisotropic prop- 
erty of the QSO emission. Both the relativistic time- 
delay effect and the anisotropy of QSO emission may 
make the apparent shape deviate from a sphere. The 
time-delay effect systematically makes the apparent size 
of the ionization front in the directions away from the 
observer smaller than the size in the directions towards 
the observer. The deviation of the shape caused by the 
time-delay effect is not significant if the QSO age is sig- 
nificantly long compared to the hydrogen recombination 
process within the ionization front and the QSO luminos- 
ity evolution is significantly slow. The time-delay effect 
does not break the rotational symmetry of the ionized 
region about the observer's line of sight to the QSO. The 
apparent expansion of the ionization front may be su- 
perluminal at the early stage of its expansion or if the 
QSO ionizing photon emission rate increases dramati- 
cally (e.g., exponentially). The anisotropic effect on the 
apparent shape depends on the detailed anisotropy of the 
emission (e.g., affected by the normal direction of the ac- 
cretion disk in the QSO) and may break the rotational 
symmetry about the observer's line of sight. 

The apparent shape of the highly ionized regions may 
be mapped by transmitted spectra of background sources 
behind or inside them or by surveys of the hyper-fine 
transition 21cm line of neutral hydrogen. The apparent 
shape could also be used to constrain the cosmological 
parameter Ha by the Alcock-Paczyhski test. 
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